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In  the  present  review  of  the  subject  special  attention  is  paid  to  the  changes  in 
turbulence  structure  brought  about  by  the  additives,  and  an  effort  is  made  to  under¬ 
stand  these  in  the  light  of  recent  advances  in  the  study  of  shear  flow  turbulence 
in  ordinary  fluids.  Also  a  simplified  theoretical  model  of  shear  flow  wall 
turbulence  is  presented,  with  the  aid  of  which  an  exploratory  study  of  the  influence 
of  some  non-Newtonian  fluid  properties  is  carried  out.  A  particularly  intriguing 
question  is  posed  by  the  remarkable  effectiveness  of  the  extremely  high-molecular- 
weight  polymers  with  which  only  a  few  parts  per  million  of  weight  of  solvent  is 
sufficient  in  some  cases  to  lead  to  drag  reductions  of  50  per  cent  or  more}  how  Is 
it  possible  that  such  minute  amounts  can  change  the  mechanical  properties  of  the 
fluid  enough  to  produce  so  dramatic  effects?  A  general  review  of  the  very  voluminous  , 
literature  on  the  subject  of  drag  reduction  has  not  been  attempted  here  (to  date 
several  hundred  publications  on  the  subject  have  appeared).  Generally,  the  emphasis 
is  put  on  the  more  recent  contributions,  particularly  those  having  a  bearing  on 
turbulence  structure,  but  some  of  the  older  materir’  of  interest  for  giving  the 
background  is  also  included. 
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1  •  Introduction 


The  reduction  of  flow  friction  drag  through  polymer  addition,  a  phenom¬ 
enon  first  subjected  to  systematic  study  by  Toms  (1948)  (with  whose 
name  it  is  commonly  associated)  together  with  Oldroyd  (1948),  but 
of  which  the  first  documented  discovery  is  due  to  Mysels  (1949),  has 
been  the  subject  of  very  extensive  research  during  the  last  decade. 
However,  despite  the  combined  efforts  of  specialists  from  many  differ¬ 
ent  fields,  no  satisfactory  explanation  of  the  phenomenon  has  as  yet 
been  found.  Drag  reduction  has  been  observed  with  a  large  number  of 
additive-fluid  combinations  (u  few  of  which  are  shown  in  figure  1), 
the  most  effective  ones  being  high-molecular-weight  linear  polymers 
(such  as  poly  ethyl  eneoxi  da )  in  water*.  There  are  clear  experimental 
indications  that  drag  reduction  is  associated  with  a  modification  of 
the  turbulence  struccurej  and  th.t3  is  most  remarkable  in  view  of  the 
permanence  and  statistical  stability  exhibited  by  ordinary  Newtonian 
shear  flow  turbulence.  Therefore,  the  phenomenon  of  drag  reduction 
becomes  of  interest  not  only  because  of  the  promising  technical  appli¬ 
cations  it  offors  but  Juit  as  much  for  fundamental  fluid  mech.-mic.il 
reasons.  Through  the  unravelling  of  the  mechanism  underlying  tho  phe¬ 
nomenon  one  should  not  only  improve  tho  prospects  of  finding  more 
efficient  and  perhaps  completely  novel  methods  of  turbulence  control, 
but  should  also  hope  to  gain  a  substantially  deepened  understanding 
of  the  fundamental  nature  of  turbulence  as  well. 

In  the  present  review  of  the  subject  special  attention  is  paid  to  the 
changes  in  turbulence  structure  brought  about  by  tho  additives,  and ' 
an  effort  is  nude  to  understand  these  ir»  the  light  of  recent  advances 
in  tho  study  of  shear  flow  turbulence  in  ordinary  fluids.  Also,  a 
simplified  theoretical  model  of  shear  flow  wall  turbulence  i3  presented, 


with  tho  aid  of  which  an  exploratory  study  of  the  influence  of  some 
non-Newtonian  fluid  properties  is  carried  out.  A  particularly  intrigu¬ 
ing  question  in  posed  by  tho  remarkable  effectiveness  of  the  extremely 
high-molecular-weight  polymers  with  which  only  a  few  parts  per  million 
of  weight  of  solvent  is  sufficient  in  some  oases  to  lead  to  drag  re¬ 
ductions  of  50  per  cent  or  more;  hew  is  it  possible  that  such  minute 
amounts  can  change  the  mechanical  properties  of  the  fluid  enough  to 
produce  so  dramatic  effects?  A  general  review  of  the  very  voluminous 
literature  on  the  subject  of  drag  reduction  has  not  been  attempted 
here  (to  date  several  hundred  publications  cn  the  subject  have  appeared). 
Generally,  tho  emphasis  Is  put  on  the  more  recent  contributions,  partic¬ 
ularly  those  having  a  bearing  on  turbulence  structure,  but  some  of  the 
older  material  of  interast  for  giving  the  background  is  also  included. 


2,  Data  on  drop;  reduction  and  moon  velocity  profiles 


Hoot  of  the  data  on  drag  reduction  published  to  date  have  been 
obtained  in  pipe  flows  following  Tons  ( 19**8)  original  procedure. 
However,  results  for  flat-plate  boundary  layers  such  as  those  by 
Kowalski  (1967)  and  by  Latto  cl  Shen  (1970)  shew  the  sane  qualitative 
and  quantitative  behaviour  as  for  pips  flow  as  do  also  those  ob¬ 
tained  in  the  rotating  coaxial  cylinder  experiments  of  Merrill  et 
al.  (i960).  This  is  perhaps  not  surprising  in  view  of  the  similar 
behaviour  of  all  turbulent  wall  shear  flows  in  Newtonian  fluids. 

We  shall  therefore  here  consider  only  the  pipe  flow  data. 

Systematic  studios  of  pipe  flow  (e.g.,  Mctzner  &  Park.  1964; 

Savins,  1964;  Elata  &  Tirosh,  1965;  Tabula,  1965;  Goren  &  Norbury, 
1967;  Seyer  &  Metzner,  196?;  Virk  et  al.,  1967  and  Virk  et  al., 

1970)  liavc  identified  several  drag  reduction  charcteristic3 
common  to  dilute  solutions  of  linear  polymers  such  as  polyethylene- 
oxide.  (Other  drag  reducing  agents  such  as  complex  soaps  with 
micellar  molecular  structures  produce  quite  different  behaviour 
with  onset  of  drag  reduction  occurring  at  very  low  shear  rate3, 
but  these  will  not  be  discussed  here. )  At  low  flow  rates  such 
fluids  obey  the  laminar  friction  law  with  a  slight  increase  in 
drag  due  to  the  increased  viscosity  caused  by  the  polymer.  At 
higher  velocities  the  flow  will  go  through  transition  at  a  Reynolds 
number  which  -  at  least  if  the  pipe  inlet  disturbances  are  mod¬ 
erately  high  -  will  not  be  appreciably  different  from  the  New¬ 
tonian  value.  Drag  reduction  sets  in  at  a  well  defined  flow  rate 
and  abovo  this  the  frictional  resistance  depends  on  the  polymcr- 
solvent  •system,  molecular  weight  and  concentration.  With  increase 
in  flow  rate  a  mximum  drag  reduction  independent  of  polymeric 


molecular  parameters  m ay  be  reached.  For  still  higher  flow  rates 
the  drag  reduction  effect  oould  eventually  diminish  due  to  polymer 
degradation.  Pipe  flow  drag  data  typical  of  dilute  poly  ethyl  eneoxide 
solutions  are  shown  In  figure  2  taken  from  Paterson  &  Abernathy 
(1970).  The  onset  of  drag  reduction  has  been  shown  by  Virk  et  al. 
(1967)  to  bo  approximately  Independent  of  pipe  diameter  and  for 
a  given  solvent-polymer  combination  to  depend  only  on  the  wall 
shear  stress.  Furthermore,  by  comparing  results  for  polymers  of 
different  molecular  weight,  they  were  ablo  to  show  that  the  onset 
data  correlated  with  tho  formula 

DMuJ/v  -  0.015  i  0.005  (1) 

v/here  is  the  effective  diameter  of  the  randotnly  coiled  macro- 
molecule  at  rest  and  u*  the  value  of  the  shear  velocity  ut  -  tjr^/p 
at  onset.  This  result  is  quite  surprising  because  It  implies  that 
drag  reduction  would  start  when  the  diameter  of  the  molecule  Is 
about  10“^  times  the  size  of  the  smallest  dynamically  significant 
eddy  in  the  flow.  Simplified  theoretical  models  for  the  macromole¬ 
cule  (PoterUn,  1966)  show  tliat  a  large  deformation  of  the  mole¬ 
cule  -  possibly  leading  to  a  complete  extension  of  the  random 
coll  -  would  take  place  when  the  shortest  time  scale  of  motion 
would  fell  below  twice  the  relaxation  time  of  the  molecule. 

However,  the  onset  data  correlate  poorly  on  a  time  hypothesis. 

Virk  (1971a)  has  also  found  that  tho  onset  of  drag  reduction  in 
rough  pipes  occurs  at  the  same  shear  stress  as  in  smooth. 

The  existence  of  a  maximum  drag  reduction  asymptote  has  by  now 

been  well  documented  for  most  known  solvent-polymer  systems. 

Figure  y,  reproduced  from  Virk  et  al.  (1970),  shows  data  presented 

1/2 

in  Prandtl's  coordinates  1/yf  versus  Re  f '  .  It  is  convenient  to 
present  vor.ults  for  pipe  flow  in  this  manner,  bocauno  in  these 


coordinates  both  the  laminar  drag  law  and  the  Prandtl  -  von  Karajan 
law  for  Newtonian  turbulent  pipe  flow  ,  v 

1//P  -  4.0  log10(Ref1/2)  -  0.4  (2) 

appear  as  straight  lines  in  a  semilogarithmie  diagram.  The  drag 
asymptote  was  also  shown  to  yield  an  approximately  straight  line 
of  the  form  (Virk  et  al,,  1970) 

! 

\/Jt  -  19.0  log,0(Re  f1/2)  -  32.4  (3) 

In  the  regime  below  the  maximum  drag  asymptote  the  friction  faotor 
varies  both  with  polymer  properties  and  flow  variables. 

Several  attempts  have  been  made  to  correlate  friction  data  in 
this  polymeric  regime  below  the  maximum  drag  asymptote  and  we 
have  singled  out  the  recent  ones  of  Seyer  &  Netxner  (1969b)  and 
of  Virk  (19/lb)  as  being  representative  of  two  different  approaches, 
one  (Seyer  &  Metzner’ s)  based  on  measured  viscoelastic  properties 
and  the  other  (Virk' a)  working  direotly  with  molecular  parameters. 
These  are  best  understood  in  light  of  +he  behaviour  of  the  mean 
velocity  profiles.  Virk  (1971b)  has  compared  published  experi¬ 
mental  data  on  mean  velocities,  and  in  figure  4  we  have  reproduced 
these  (figure  2  of  Virk,  1971b)  showing  the  measured  /elo cities 
plotted  in  the  usual  wall  variables  U+  =>  U/ut,  yu^/v.  Prom  those 
one  can  distinguish  three  different  regions  which  can  bo  repre¬ 
sented  approximately  by 


(I) 

U+  -  y+  i  0  <  jr* 

+ 

<yv 

(4) 

(II) 

U+  -  Agin  y++  Dg 

•I*  +  + 

»  y/y  <yc 

(5) 

(III) 

U+  «  A.^ln  y+  +  By 

1  y  e  <  y+  * 

(6) 

respectively,  whore  R+  is  the  value  nf  y+  correspondin';,  to  the 
.  .  .  ftn  is  th.O  USU!il  V'iSCOUS  SUblO„  0» 
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cxtendln;*  to  about  y+  -  10.  The  second  region,  called  the  "clastic" 

region,  :1s  best  i»:pi>r.:.r.tc:l  by  the  values  Ag  -  11.7  and  llg  --  -  17.0. 

In  the  third  .’onion  t/  shows  the  name  logarithmic  variation  with  y' 

an  in  the  Hcwtonlon  case  but  the  cui've  in  displaced  upwards.  Thus, 

A..  -  K.b,  but  11.  Is  a  function  of  molecular  and  flow  parameter. 

J>  J 

The  mean  velocity  profile  data  hence  indicate  the  existence  of  a 
Newtonian  "plus"  3n  the  center  of  the  pipo  with  a  mixing  length 
apparently  unaffected  by  the  polymers.  At  the  condition  of  maximum 
&c?z  reduction  the  clastic  region  II  is  assumed  to  extend  all  the 


way  to  the  center  of  the  pipe,  and  the  Newtonian  region  III  then 
disappears.  Vii'k  (1971b)  dctci mired  the  coefficients  in  the  region  II 
so  as  to  make  the  profile  give*  a  drag  variation  consistent  with  the 
maximum  dreg  asymptote  "odor  this  assumption.  Using  (’0  -  (6)  one 


can  integrate  the  total  volumetric  flov;  rate  and  hence  obtain  a  re¬ 
lation  between  friction  force  and  flow  rate. 


Seyer  and  Mctzncr  (1969b)  assumed  the  viscous  wall  layer  to  extend 
all  the  way  to  the  Newtoni  an  region  without  an  intc  'mediate  region.  . 
Thus,  setting 

U+  *  Alny++B(et)  (7) 

where  d+  -  Ou2  /v  ond  0  is  the  relaxation  time  of  the  fluid  they 

T 

wore  able  to  correluto  data  for  different  polymer  solutions.  The 
constant  A  is  the  same  as  in  the  Newtonian  case,  (Ueyer  and  Metznor 
took  A  «  2J\ 6.)  The  nondimcnsionol  relaxation  time  0h  is  the  Deborah 
number  formed  by  the  ratio  of  fluid  relaxation  time  to  natural  vis¬ 
cous  time  3calc  of  the  turbulence.  The  curve  of  the  intercept  func¬ 
tion  B(0+)  obtained  in  thi3  manner  is  reproduced  in  figure  5. 

Their  data  indicate  an  asymptotic  value  of  B  of  about  JO  for  large 
relaxation  tinea.  Virk  (1971b)  made 


a  comparison  of  the  Seyer  -  Metzner  model  with  his  newer  one  and 
found  that  the  latter  is  superior  in  representing  measure  profiles 
at  conditions  of  large  drag  reduction.  The  Seyer  -  Metzner  cor-  ; 
relation  is  also  difficult  to  use  for  highly  dilute  solutions, 
for  which  no  reliable  methods  of  determining  the  relaxation  time 
6  are  as  yet  available  (see  further  the  following  seotion). 

Virk  (1971b)  tries  to  relate  the  intercept  constant  in  the 
following  manner 

Bj  -  Bn  +  6  In  (R’/ft4*) 

where  Bn  is  tho  value  for  a  Newtonian  flow,  B  -  f>,5.  The  factor 
ln(R  /R  ),  where  R  is  the  value  of  R+  at  onset,  is  assumed  to 
represent  an  "excitation"  factor  which  depends  both  on  flow'  para¬ 
meters  and  such  polymetrie  parameters  that  may  be  associated  with 
the  excitation  of  the  macromolecule.  The  function  6,  wlJLch  Virk 
tenues  a  "olopo  increment  factor"  (since  it  gives  the  fractional 
Increase  due  to  tho  polymer  in  the  slope  of  tho  curve  of 
versus  log (Re  Jf)  )  is  assumed  to  depend  3olbly  on  molecular 
parameters.  By  examining  a  large  body  of  data  Virk  (1971b)  estab¬ 
lished  that  for  each  polymer- solvent  system  6  is  approximately 
pronominal  to  c  ,  vhero  c  is  the  weight  fraction  of  polymer  j 
addit  ,  Sotting  i 

I 

j 

6  «  k^(Ko/U)  1/Sf  (Q) 

whore  A  is  Avogadros  number,  c  the  weight  fraction  and  M  tho  j 

molecular  weight,  respectively,  of  the  polymer  (so  that  the  j 

j 

quantity  under  square  root  represents  the  number  of  macromoleculea  1 

per  unit  weight  of  solution)  Virk  find3  that  for  all  polymer  j 

r 

species  with  tho  exception  of  one  (guar  gum)  the  function  If  i‘ 

| 

(turned  "Intrinsic  slope  in  cement  factor")  cun  bo  fairly  uol  1 
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rcprosented  by  a  single  curve  with  the  number  of  macro  molecular 
baokbone  chain  links,  N,  as  a  parameter.  His  results  are  repro¬ 
duced  in  figure  6  (Virk' s  figure  5) .  Furthermore,  in  a  logarithmic 
diagram  all  the  data  with  tho  exception  of  that  for  guar  gum  were 
found  to  cluster  around  a  single  straight  (line  of  slope  jj  indi¬ 
cating  a  simple  power  law  of  the  form 


ir-KN5''2 


v/hero  K  is  a  constant  found  to  have  the  approximate  value 

-14 

K  »  (2.3  ±  0,8)  x  10  .  From  a  theoretical  reasoning  based  on 

the  estimate  of  the  turbulent  strain  energy  stored  in  the  flexible 
macromolceule,  Virk  comes  to  the  conciliation  that  Tf  should  be 

3 

proportional  to  b  ,  where  b  is  the  bond  length  cf  the  nomomcr 
making  up  the  macromolecule.  Guar  gum  with  u  bond  length  of 
b  ~  20  A,  as  compared  to  the  other  polymers  having  b  ~  4  A, 

3 

should  therefore  bo  about  (20/4)  as  effective  as  the  others, 
as  indued  roughly  indicated  by  tho  data  ior  guar  gum  which  give 
Kc«2fl0x1<fU. 


3«  Rheological  properties  , "f 

,  l 

Ono  of  the  most  important  and  intricate  questions  in  drag  reduction 
research  has  been  to  what  specific  bulk  rheological  property  or 
properties  of  the  fluid  the  effect  might  be  ascribed  to.  It  is  1 

usually  stated  (see,  e.g.,  Metzner  &  Park,  1964)  that  drag  reducing  i 

i' 

fluids  show  viscoelastic  behaviour,  and  the  suggestion  that  visco¬ 
elasticity  is  the  cause  of  drag  reduction  has  been  almost  uni- 
verjiolly  accepted.  Indeed,  in  larger  concentrations  the  polymers 
effective  in  reducing  the  drag  produce  normal  stress  differences 
in  pure  shear  of  moderate  magnitude,  an  effect  usually  associated  jj 
with  viscoelasticity,  Direot  demonstration  of  reversible  elastic 

jj 

effects  has  also  been  given  by  Darby  (1970),  who  with  the  aid  of  f! 
a  specially  designed  transient  rheometer  was  able  to  measure  ' 

elastic  rebound  for  fluids  with  polymer  concentrations  of  100  i 

to  1300  parts  per  million.  However,  such  rheometric  data  are  usu¬ 
ally  taken  at  low  shear  rates  for  which  the  normal  stresses  are  !'1 

♦ 

i 

expected  to  be  much  lower  than  at  the  high  shear  velocities  of  ;j 

1 

interest  in  drag  reduction.  Normal  stress  differences  for  drag  I 

!  I 

reducing  fluids  are  therefore  usually  determined  from  the  measure-  ! ; 
ment  of  the  thrust  of  a  capillary  Jet  (see  Metzner  &  Park,  1964).  i 

J  u; 

a  method  which  does  not  suffer  from  this  shortcoming.  This  method  i  f  ’ 

•  r 

wan  the  ono  used  by  Oadd  (1966a,b)  and  by  Oliver  (1966),  the  j ! 

latter  also  reporting  results  for  fairly  low  concentration  levels 

f  I 

at  which  the  method  becomes  unreliable,  however.  With  the  aid.  of  l! 

!  j 

a  specially  designed  concentric  tube  rheometer,  Oadd  (1966b),  was  f’ 
e\so  able  to  detect  a  non-zero  second  normal  stress  difference 
^0V  p.p.n, )  poly  ethyl  cncoxidc  solution 

but  did  not  find  any  appreciable  second  normal  stress  difference 


acrylamide. 


It  has  also  been  found  that  the  viseoelabtio  properties  exhibited 
depend  strongly  on  aging  and  degradation  of  the  polymer  solution, 
especially  in  the  cr.se  of  polyethyleneoxide.  Brenr.cn  end  Oadd 
(1967)  reported  a  much  diminished  total  hoad  loss  with  age  of 
solution  for  a  pitot  tube  in  a  stream  of  polyethyleneoxide  solu¬ 
tion  indicating  a  diminished  first  normal  stress.  They  also  stated 
that  the  socond  normal  stress  for  a  70  p.p.ra.  solution  vanished 
with  aging.  Several  later  investigatiors  have  noted  the  anomalous 
flow  behaviour  of  fresh  polyehtyleneoxide  solution}  thus 
Kalashnikov  and  Kudin  (1970)  found  that  the  Strouhal  number  for 
tho  Karman  vortex  sheet  behind  a  circular  cylinder  was  substanti¬ 
ally  reduced  in  a  flow  of  fresh  solution,  whereae  the  same  solu¬ 
tion  passed  through  a  eentrifugal  pump  gave  a  small  increase  in 
frequency  a3  compared  to  water,  as  did  a  solution  of  guar  gum. 
Despite  such  marked  degradation  effects  a  degraded  solution  was 
still  able  to  produce  a  'ubstnntial  drag  reduction  (although 

• 

somewhat  diminished  as  compared  to  that  of  a  fresh  solution). 

Such  behaviour  has  led  Barenblatt  et  al.  (1965)  to  propose  that 
the  viscoelastic  behaviour  of  a  polyehtyleneoxide  solution  is 
a  result  of  tho  formation  of  very  large  aggregates  of  macro- 
molocules  w'vich  ir  part  disappear  with  aging.  The  larger  aggre¬ 
gates  will  also  bo  broken  up  by  the  shear  field  in  tho  flow 
apparatus,  00  that  lnrGo-ocnle  motion  such  as  vortex  shedding 
would  bo  affcotod  by  degradation  but  not  the  drag  reduction. 

Although  standard  rhcometric  measurements  for  very  dilute  polymer 
solutions  (<  20  p.p.m. )  usually  Give  no  clear  indication  of  depar¬ 
ture  from  a  purely  Newtonian  liquid  (Bryson  et  al.,  1971 ),  there 


.11. 


is  other  evidence  that  such  .fluids  can  show  substantial 
non-Newtonian  behaviour.  Metzner  and  Metzner  (1970)  measured  the 
stress  levels  developed  In  dilute  solutions  of  polyacrylamide 
Issuing  from  small  orifices  at  high  velocities.  In  this  flow  field 
the  deformation  is  approximately  that  of  pure  rapid  elongation. 

The  ratio  of  olongational  to  shear  viscosity  (which  Is  equal  to 
three  in  a  Newtonian  fluid)  determined  in  this  way  was  found  to 
be. exceedingly  high;  a  conservative  estimate  of  the  deformation 
rate  level  yielded  values  as  high  as  29,000.  Possibly,  the  ob¬ 
served  tendency  of  polymer  solutions  to  form  long  threads  may  be 
related  to  such  high  values  of  bulk  viscosity  (Lunacy,  1969). 

In  order  to  use  laboratory  rheological  data  in  tho  equations  of 
motion  for  a  continuum  fluid  one  needs  to  incorporate  them  in  a 
constitutive  model  for  the  fluid.  A  large  number  of  such  models 
havo  been  proposed  and  studied  in  the  rheological  literature. 
Spriggs  et  el.  (1966)  mrdc  en  appraisal  of  twenty  different  models 
for  viscoelastic  materials,  end  they  concluded  that  in  order  to, 
reproduce  all  essential  rheological  properties  a  model  with  at 
least  four  different  constants  was  needed.  In  view  of  the  uncer¬ 
tainty  as  to  how  viscoelastic  effects  come  into  play  In  the 
phenomenon  of  drag  reduction,  however,  it  would  seem  most  approp¬ 
riate  to  select  for  exploratory  studies  the  simplest  model  incor¬ 
porating  viscoolasticity.  Such  a  model  is  tho  convectod  Maxwell 
one  proposed  by  White  and  Metzner  0965),  which  in  Cartesian 
coordinates  may  bo  written 

Tu +  0  wl  ■  is  -  ■  ('0) 
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where  t 


io  the  deviatoric  shear  stress,  0  a  relaxation  time. 


6/ fit  the  Oldroyd  time  derivative 
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and  and  d^  have  their  usual  meaning  of  shear  viscosity 
and  deformation  rate  tensor 
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respectively*  The  material  coefficients  6  and  p  generally 
depend  on  the  deformation  rate*  Tho  constitutive  relation  (10)  is 
a  simplified  version  of  Oldroyd' s  (1958)  original  one  containing 
eight  material  constants,  for  exploratory  purposes  it  is  a  great 
advantage  to  have  a  model  with  only  one  posltivo  constant,  namely 
0,  to  describe  departure  from  a  Newtonian  fluid.  It  has  also  some 
foundation  in  physical  reality  as  Luraley  (1971)  has  shown  that  a 
simple  model  for  a  dilute  solution  of  polymer  molecules  employing 
a  dumbell  model  for  the  molecules  leads  to  a  constitutive  relation 
of  the  Oldroyd  type  (albeit  with  three  instead  of  two  material- 
coefficients).  It  is  also  the  relation  usually  employed  to  con* 
vert  data  on  measured  normal  stress  differences  to  relaxation 
times j  application  of  (10)  to  simple  steady  3hoar  flow, 

U  "(UjfXgJAO),  gives 

T  j  j  *  T  22  u  *  00 


whoro  T  “  SUj/ftXg.  In  this  model  the  fluid  relaxation  time  0 
is  thus  proportional  to  the  first  normal  stress  difference*  On 
tho  basis  of  (10)  ono  can  also  demonstrate  the  possibility  of 
high  resistance  duo  to  stretching  as  observed  in  the  experiments 
by  Mcrtzncr  and  Kotanov  (1970)  quotod  above* 


In  view  of  the  great  difficulties  In  measuring  departures  from 
Newtonian  behaviour  of  very  dilute  drag  reducing  fluids,  it  would 
be  natural  to  turn  to  theory  to  see  what  guidelines  it  may  offer. 

In  view  of  tho  large  size  of  the  macromoleeules  in  comparison 
with  tho  solvent  molecules  involved  a  good  mechnnistic  model  for 
a  drag  reducing  fluid  is  to  recard  it  as  a  suspension  of  small 
particles  in  a  continuum  fluid.  Tho  rheology  of  suspensions  is 
a  field  of  research  which  has  undergone  a  vcr,y  extensive  develop¬ 
ment  since  its  inception  by  Albert  Kinsteln  in  lyO 'j  with  Ms 
famous  analysts  of  the  increase  in  fluid  viscosity  duo  to  the 
presence  of  rigid  spheres.  Good  review  articles  on  tho  subject 
con  be  found  in  Vole.  I  and  IV  of  the  series  "Rheology",  edited 
by  fed  rich  ( 1956,  1966).  In  theoretical  suspension  mechanics  one 
can  distinguish  the  following  three  separate  stages  in  tho  ana¬ 
lysis: 

i)  calculation  of  tho  influence  of  the  local  flow  field  on  the 
individual  particle) 

ii)  determination  of  tho  effect  of  a  statistical  ensemble  of 
particles  on  the  bulk  rheological  properties  for  the  suspension)  and 

iii)  analysis  of  the  bulk  motion  of  the  suspension. 

In  i)  fluid  inertia  can  bo  ignored  because  of  the  smell  size  of 
the  particle,  but  the  calculation  may  nevertheless  become  very 
complicated  if  the  purticlo  undergoes  large  deformations  or  is 
subjected  to  hydrodynamic  interference  from  neighboring  particles, 
situations  that  do  arise  for  the  fluid  of  interest  hers.  Also, 
tho  effect  of  Brownian  motion  may  need  to  bo  incorporated  in 
determining  the  shape  ond  averago  orientation  of  tho  particles, 

Tho  greatest  complication  is  perhaps  caused  by  the  coupling 
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fa  etween  stages  ill)  and  i)  j  the  bulk  motion  of  tho  suspension  Is 
Influenced  by  tho  chape  and  orientation  of  the  particles,  which 
In  tum  depend  on  the  local  flow  field  created.  The  problem  of 
determining  the  shape  of  the  macronolecule  can  also  become  quite 
difficult  depending  on  the  degree  of  sophisti cation  with  which  the 
molecule  is  modelled.  The  11m  ar  co-call cd  bead-  and- spring  models 
employed  by  Rouse  0952)  and  by  Zimm  (1956)  have  been  employed  to 
calculate  the  relaxation  time  of  a  dilute  suspension  of  coiled 
maorontolccules  subjected  to  a  small  oscillatory  shear  (James  & 
Acosta,  1070) .  i’ctorlJn  (i *>>C)  employing  a  linear  dumbo'll  model 
(l,c,,  two  beads  and  one  spring)  demonstrated  tnat  tho  molecule 
can  become  infinitely  extended  Jn  an  irrotational  straining  field 
if  the  strain  rate  exceeds  the  inverse  of  twice  the  relaxation 
time  of  the  molecule.  In  a  pure  shear,  on  the  othnr  hand,  it  will 
never  become  fully  extended. 


Tho  infinite  extension  is  of  course  a  consequence  of  the  t-scumpelon 
of  spring  linearity.  V/Jth  the  aid  of  a  more  reaJ.isl.tr.  nonlinear ^ 
thread  model,  Which  (1°72)  has  shown  recently  that  there  in  a 
critical  strain  r. Wo  a’ncvc  which  the  molecule  becomes  compxotcy 
(at  least  ’./*•)  extruded}  below  this  tho  molecule  is  not  extended, 
at  all.  One  would  therefore  expect  that  in  a  turbulent  field  of 
su.fi.te.'.-'  at  intensity  a  groat  proportion  of  the  macro.-, olcculcn  *3 
uncsilc-1.  On  the  of  Rti.r-Un'r.  (1'j'9'>)  result  La  rUy  (V.)C')) 

tinted  that  7C7>  of  all  molecules  should  be  extended  in  a  tur¬ 
bulent  flow  if  tho  root  mo  in- square  strain  rate  oxecodr.  tha  oriti- 


cal  value  of  the  Inverse  of  twice  the  relaxation  t T.; so. 


;;uoh  simple  rheological  theories  point  to  the  im¬ 
portant  conclusion  that  the  molecular  conformation  that  may  be 
vay  common  in  a  turbulent  flow  might  be  an  extended  thread 
rather  than  a  random  coil,  go  that  the  rheological  properties 
of  the  polymer  solution  in  a  turbulent  flow  therefore  might  bo 
quite  different  from  what  they  are  in  a  laminar  flow.  As  a  theory 
realistically  modelling  coiling  -  uncoiling  m»cromolecules  in  an 
unsteady  flow  field  is  likely  to  lead  to  un  surmount  able  mathe¬ 
matical  difficulties,  it  would  therefore  be  more  useful  for  an 
exploratory  study  to  search  for  a  model  which  gives  the  maximum 
rheological  effects  of  extended  partieles. 


Such  a  model  can  in  fact  be  found  from  the  extensive  recent  studies 
of  the  rheological  properties  of  suspensions  of  rigid  elongated 
particle3  carried  out  by  Batchelor  (197<>,b,  ’971 )  and  his  group 
at  Cambridge  University.  The  main  concern  in  this  group  has  been 
with  highly  elongated  rigid  particles.  The  problem  of  determining 
the  partlclo  motion  that  simplifies  to  that  of  determining  the  . 
velocity  and  orientation  giving  zero  not  force  and  moment  on  the 
particle.  (Inertia  of  the  particle  is  neglected.)  This  was  a 
problem  considered  by  Jeffrey  0922),  who  determined  the  oi’bit 
of  spheroids  in  a  shearing  motion  and  in  an  trrotational  straining 
field.  Leal  &  Hinch  (1971)  and  lilnch  &  I.cal  (1972)  have  extended 
this  analysis  to  incorporate  Brownian  motion  and  have  shown  that 
the  orbit3  of  highly  elongated  particles  in  shear  flow  tend  to 
become  ouch  as  to  make  them  spend  moot  of  their  time  with  thoir 
longest  axes  aligned  with  the  direction  of  flow.  Jeffrey  (1922) 
al.'-o  demon?’. tratM  that  in  rn  irrotational  strMntng  motion  all 


the  pfrtlcj.es  would  end  up  in  a  direction  with  thoir  longest  axes 
a]  v-nrd  with  tjv  flow.  limee,  there  An  a  strong  tendency  for 


<•1  :>.i' 


to  1.1.}';:  up  in 


flow  dlr'^-tlon.  If 
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ticles  are  very  slender,  their  contribution  to  the  stress  in  pure  ’ 
shear  will  be  quite  small  so  that  the  suspension  , 

will  behave  nearly  like  a  Newtonian  fluid.  On  the  other  hand,  any 
3 training  perturbation  in  the  mean  flow  direction  would  be  strongly 
resisted.  We  can  use  such  general  results  together  with  that  for  the 
particle  stress  contribution  due  to  pure  straining  motion  calculated 
by  Batchelor  (1970  to  determine  the  perturbation  stress  t,,  arising 
from  on  unsteady  twodimensioned  perturbation  with  velocity  components 
u,  -  u  and  u2  «  v  of  a  parallel  shear  flow,  a  model  problem  that,  as 
will  be  argued  below,  has  some  bearing  on  the  dynamics  of  a  turbulent 
shear  flow. 

Batchelor  showo  that  the  contribution  At,,  to  the  stress  in  pure 
straining  motion  due  to  the  presence  of  slender  particles,  all  of 
length  1  end  diameter  d,  is  given  approximately  by 

-  BU1 

4Tu  '  •  <'2> 

where 

5  -  |  (  r2/ln(n/| )  «  2  nl?/ln(rr/« )  ,  (1?) 

ft  is  the  volume  fraction  of  particles,  r  =  l/d  their  aspect  ratio, 
and  n  their  number  density,  for  a  perturbation  in  the  mean  flow  direc¬ 
tion,  the  perturbation  stress  would  therefore  be  obtained  by  the  above 
formula  with  U,  replaced  by  u.  However,  the  perturbation  velocity 
component  v  normal  to  the  flow  would  in  addition  cause  a  misalignment 
of  the  fluid  particles  with  the  (small)  angle 

whoro  the  symbol  Dt  means  that  the  integration  should  be  carried  out 
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following  tho  particle.  Since  the  partloles  are  now  tilted  with  respect 
to  tho  mean  flow  they  will  experience  a  velocity  gradient  along  their 
length  proportional  to  c(dU/dy)j  henoe  the  total  first-order  contri¬ 
bution  to  the  stress  will  be 


t 


This  result  has  been  derived  by  Hinch ,  who  also  gave  estimates  of  the 
error  arising  because  not  all  particles  are  aligned.  We  will  in  the 
following  refer  to  (15)  as  the  Batchelor  -  Hinch  model.  Adding  the 
contribution  due  to  the  solvent  itself,  which  is  obtained  from  the 
ordinary  Newtonian  formula,  one  will  thus  find  the  total  perturbation 
stress  due  to  the  perturbation  velocity  field.  This  result  will  be 
U3cd  below  wi  n  a  simplified  model  for  a  turbulent  shear  flow  exploring 
lt3  large-  and  small-scale  stability  properties,  Batchelor  ( 1971 ) 
derived  the  formula  (12)  for  rigid  rods  under  the  restriction  that 
the  suspension  is  dilute  in  the  volumetric  sense  so  that  the  average 
distance  between  the  rods  is  large  compared  to  their  diameter  but  that 
the  rod  length  is  much  larger  than  the  average  distance  between  the 
rods  as  illustrated  in  figure  7»  Under  such  conditions,  wliich  would  be 
satisfied  for  the  dilute  polymer  solutions  of  present  interest  if  tho 
macromolecular  coils  were  extended,  the  coefficient  B  con  be  very  largo, 
indeed,  .leading  to  very  largo  ratios  of  extensicnal  viscosity  to  shear 
viscosity  of  the  order  of  10'^  or  more,  values  that  order- of- magnitude  - 
wise  are  consistent  with  the  exit eri mental  results  reported  in  the  paper 
of  Metzner  and  Motzner  ( 1970)  quoted  above.  Where  the  model  obviously 
fails  is  in  the  assumption  of  rigid  rods,  A  fluid  filled  with  threads 
# 

K.J.  Hinch,  private  communication. 
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from  uncoiled  mccromolecules  would  be  able  to  resist  extension  but 
not  compression.  A  more  reallstlo  model  for  such  a  fluid  would  have 
to  be  nonlinear  and  quite  complicated,  however,  and  in  view  of  the 
fact  that  some  suspensions  of  elongated  inelastic  particles  have 
shorn  drag  reduction  (Bobkovdcz  &  Gauvinj  1965),  an  investigation 
based  on  the  simple  Batchelor  -  Hinch  model  (15)  should  be  of  some 
value  for  an  exploratory  study* 


it»  Stability  and  transition 


Tn  the  mechanics  of  turbulence  the  phenomena  of  instability  and 
transition  of  a  laminar  flow  are  believed  to  be  of  fundamental  impor. 
tance.  Therefore,  many  theoretical  and  experimental  investigations 
have  been  aimed  at  a  study  of  non-Newtonian  effects  on  instability 
and  transition. 

The  theory  of  hydrodynamic  stability  is  concerned  with  whether  small. 

disturbances  in  a  steady  larlnar  flow  field  will  grow  or  decay.  We 

shall  here  discuss  only  two-dimensional  infinitesimal  disturbances 
of  a  piano  parallel  shear  flow. 

Ui  "  (u(*2)A°)  (16) 

(Xj  »  x,  "  y,  Xj  -  z) 

For  such  disturbances  it  is  a  straight- forward,  if  often  somewhat 
tedious,  process  to  derive  from  the  Unearned  momentum  and  the 
chosen  constitutive  relation  a  single  differential  equation  for  the 
stream  function  of  the  perturbation.  Setting  in  the  usual  manner  for 
the  perturbation  velocity  components 

u  -  (Dc?)  elQr(x"cb)  (17) 

v  -  -  i(v<pelc^x-ct) 

where  a  and  c  arc  the  wave  number  and  phase  velocity,  respectively, 
of  the  »;avo-llko  disturbance}  n  .  d/dy;  and  the  real  part  of 
t,o(y)exp  [ior  (x-fct)  j  is  the  perturbation  stream  function,  one  obtains 
for  the  two  constitutive  models  selected  the  following  modified 
Orr  -  Sommorfnld  equations t 


a)  Convectcd  Maxwell  model: 


(U-c)(D2-(y2)9-9D2U  -  [(D2-  a2)?\  +  \  hnD2(,l 

L  n«0  J 


(18) 


where  Re  is  the  Reynolds  number,  and 
S  »  1  +  iaQ(U-c) 

bO  “  “  or2(D2S)T  -  2*2(D3)2(1  +if2) 

-  4(DS)(iy5S)3"1-  jXD^jV1  v  4(DS)V2 

-  6(DS)2(D23)T^1 

bj  »  -  2c?(DS)T  +  2n?sr  +  4(ds)(d2s)t2 
-  4(ds)^ts”^ 
b2  »  3D2ST-i-2(DS)2T 


b..  «  2D  ST 
5 


This  equation  was  given  by  Porteoua  and  Dem  (197I),  who  also  demon¬ 
strated  that  for  small  fluid  elasticity,  i.e.,  small  6,  only  the 
underlined  term  In  the  above  expressions  for  should  bo  retained. 
The  equation  then  becomes  equivalent  to  the  one  for  a  second-order 
fluid  employed  by  Chun  &  Schwartz  (1968)  and  others,  an  equation 
also  arrived  at  by  Walters  (1962)  by  considering  a  fluid  of  short 
memory, 

b)  Batchelor  -  Hinch  model: 

(U-o)  (D8-  „8  )  7  -  ?  \h  .  -jj;-  [(D8..  „2)8  9 
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All  quantities  appearing  in  (18)  and  (19)  have  been  nondimensionalizea 
in  the  usual  minor  for  boundary  layers  employing  shear  layer  thickness 
6  and  free  stream  velocity  as  reference  length  and  velocity,  respec¬ 
tively,  end  6/U1  as  reference  time  scale.  The  Batchelor  parameter  B 
is  nondimcnsional  as  such  and  is  independent  of  flow  parameters. 

Moot  stability  results  available  to  dato  have  been  based  on  the  second 
order  fluid  constitutive  model.  For  plane  Poiseulle  flow  Chan  Man  Fong 
and  Walters  (1965)  found  that  a  small  elasticity  tended  to  lower  the 
critical  Reynolds  number.  This  destabilizing  effect  way  confirmed  by 
Chun  &  Schwartz  (1968)  through  direct  numerical  calculations  of  disturb¬ 
ance  amplification  rates,  which  for  moderate  values  of  the  elasticity 
parameter  were  found  to  increase  somewhat  o,rcr  those  for  a  Newtonian 
fluid.  Further  demonstrations  of  the  destabilization  due  to  elasticity 
has  also  been  given  by  Porteous  and  Denn  (1971)  who  employed  Doth  (18) 
and  the  second  order  fluid  model  and  found  very  little  difference 
between  the  results  for  small  values  of  the  fl’ld  elasticity.  They  also 
carried  out  an  analysis  for  finite  amplitudes.  No  results  obtained  cn 
basis  of  the  Batchelor  -  Hinch  model  (19)  for  aligned  rigid  rods  have 
as  yet  published,  but  our  own  preliminary  calculations  reported  below 
indicate  that  their  effect  should  generally  b«  stabilizing.  In  the 
literature  there  have  also  appeared  several  analysis  of  stability  of 
Couettc  and  Taylor  -  Couctto  flow  of  viscoelastic  fluid,  of  which  per¬ 
haps  the  most  ambitious  one  as  far  as  the  constitutive  model  employed 
in  that  of  Kavlsson  et  al.  (1971).  Elasticity  is  found  sometimes  to  act 
stabilizing  for  a  Taylor  -  fouettc  flow  but  it  is  difficult  to  Judge 
what  possible  implications  such  results  may  have  for  shear  xlow 


turbulence. 


Experiments  directly  aimed  at  studying  the  stability  of  drag  reducing 
fluids  have  been  carried  out  by  Berman  &  Cooper  (1972),  They  investi¬ 
gated  the  flow  response  due  to  the  periodic  motion  of  a  thin  sleeve. 
Although  the  shear  rates  were  not  high  enough  for  viscoelastic  effects 
to  be  expected,  they  found  that  the  periodic  disturbance  generated 
by  the  sleeve  persisted  longer  downstream  with  the  dilute  polymer 
solution  than  vdth  the  pure  solvent. 

The  pipe  flov/  transition  experiments  that  have  been  reported  in  the 
literature  to  date  give  rather  contradictory  information  regarding 
whether  p  drag  reducing  fluid  will  speed  up  or  delay  transition  to 
turbulence,  Virk  ot  al.  (1967)  reported  no  substantial  effect  of 
polymer  addition.  Castro  and  Squire  (1967)  and  White  &  MeEligot  P970) 
found  a  slight  increase  in  the  transition  Reynolds  number,  whereas 
Chung  &  Graebel  (1972)  obtained  a  large  increase.  Paterson  and 
Ahomathy  (1972),  on  the  other  hand,  obsei'ved  a  decrease  in  the 
transition  Reynolds  number  of  strongly  drag-reducing  solutions  as 
compared  to  the  pure  solvent,  ‘ 

Possibly,  those  strongly  differing  results  could  in  part  have  their 
explanation  in  the  different  methods  rf  detecting  transition  used. 
However,  it  is  also  well  known  from  pipe  flow  transition  measurements 
in  Newtonian  fluids  (Lindgren,  1997)  that  by  reducing  inlet  distur¬ 
bances  ono  can  prolong  tho  laminar  regime  almost  indefinitely.  There¬ 
fore,  one  can  expect  that  only  experiments  conducted  at  reasonably 
high  inlet  disturbance  levels,  such  as  those  of  Virk  et  al,  (I967) 
will  yield  reproducible  results  for  the  lower  transition  Reynolds 
number.  It  is  of  interest  that  .Paterson  &.  Abernathy-  (1972),  when 
using  a  squared-off  pipe  producing  disturbed  inlet  flow,  found  little 
effect  from  tho  polymer  on  the  transition  Reynolds  number,  in  agree¬ 
ment  with  t.ho  findings  of  Virk  et  el,  (1967) . 


5.  Turbulence  structure 


In  recent  years  considerable  advances  have  been  made  in  the  under¬ 
standing  of  the  dynamical  processes  taking  place  in  a  turbulent 
Newtonian  wall  shear  How.  Earlier  measurements  of  the  turbulent 
field  were  focussed  on  determining  space-time  correlations  of  the 
fluctuating  quantities.  From  such  long-term  averages  one  could  hope 
to  detect  any  persistent  or  recurring  larger-scale  flow  features  that 
might  be  of  dynamical  significance.  Measurements  of  correlations  3uch 
as  those  of  Fovro  ct  al.  (1957)  and  of  Willmarth  &  Tu  (1967)  depict 
the  main  structure  as  a  decaying  and  convected  pattern  which  propagates 
with  a  velocity  close  to  but  not  exactly  equal  to  the  local  flow  velo¬ 
city.  Corcos  (1964)  carried  out  a  spectral  analysis  of  the  wall  pressure 
fluctuation  correlations  obtained  by  Willmarth  &  Wooldridge  (1962)  and 
showed  that  th^  convection  velocity  had  a  uniquely  defined  value  for 
each  frequency.  This  was  interpreted  as  evidence  for  the  presence 
of  waves  in  the  turbulent  boundary  layer  by  Landahl  (1967),  who 

showed  theoretically  that  in  a  system  admitting  propagation  of  waves. 

* 

(a  wave  guide)  the  least  damped  wave  mode  should  dominate  in  the 
statistical  average,  Morrison  &  Krone, uer  (1969)  detected  such  waves 
from  measurements  of  frequency-wavc-num’oor  spectra  of  the  fluctu¬ 
ating  velocities  in  the  wall  layer  of  a  turbulent  pipe  flow  and  found 
that  the  fluctuating  velocity  field  is  dominated  by  waves  that  are 
highly  swept,  Hussain  £  Reynolds  (1972)  excited  periodic  disturbances 
in  a  two-dimensional  turbulent  channel  flow  mid  measured  their  ampli¬ 
tude  distributions  and  propagation  constants.  They  shewed  that  the 
effect  of  background  turbulence  must  be  included  in  the  equations 
modelling  the  disturbances  (Reynolds  C.  Hussain,  1972).  An  entirely 
different  approach  in  analysing  correlation  data  was  taken  by  Balccv.’ell 
&  Lumlcy  (1967),  who  decomposed  the  correlation  functions  in  a  set  of 
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orthogonal  eigenfunctions  and  demonstrated  that  the  dominating  l&xtge- 

scale  feature  in  the  wall  layer  consists  of  a  pair  of  randomly  ocur- 

ring  contrarotating  streamwlso  vortices  separated  in  the  spanvdse 

direction  of  a  distance  of  about  X*  ■  100  expressed  in  wall  vart- 

z 

ables. 

The  conceptual  two-layer  model  proposed  by  Townsend  (1956),  which 

divides  the  turbulent  boundary  layer  into  an  "active"  part  comprising 

the  region  closest  to  the  wall  and  a  "passive"  outer  one,  hits  been 

a  source  of  inspiration  for  many  of  the  later  invistigationn.  That 

the  energetics  processes  in  the  boundary  layer  are  mainly  confined 

to  the  inner  layer  near  the  j-rall  was  given  detailed  confirmation  by 

Kline  et  al.  (1967)  and  by  Kim  et  al,  (1971).  Their  visual  observations 

of  a  thick  turbulent  boundary  layer  in  water  with  the  aid  of  hydrogen 

bubble  tracers  revealed  that  most  of  the  turbulent  energy  and  Reynolds 

stress  was  produced  during  intermittently  ocurring  localized  "bursts" 

or  turbulence.  The  burst  is  a  violent  instability  process  found  to 

occur  primarily  in  the  zone  0  <  y+  <  100,  The  observations  suggested 

% 

a  local  structure  like  a  hair-pin  snaped  vortex  lifting  up  from  the 
lower  part  of  the  boundary  layer  and  being  stretched  by  the  mean  shear 
as  it  is  transported  downstream  (see  figure  8  reproduced  from  Kline 
et  al,,  1967).  The  average  transverse  distance  between  suen  events 
was  found  to  be  rbout  X*  100  expressed  in  wall  variables,  strongly 
suggesting  a  connection  to  the  statistical  picture  of  Bakewcll  and 
Lumley  (1967) .  Further  visual  observations  by  Corlno  and  Brodkey 
(1969)  revealed  two  distinct  violent  processeaj  an  "ejection"  of  low- 
speed  fluid  away  from  the  wall  and  a  "ev/cep"  of  high-speed  fluid  to¬ 


wards  the  wall. 
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Decause  of  their  strong  intemittency  the  presence  of  turbulent 
bursts  cannot  normally  be  detected  by  means  of  ordinary  statistical 
methods,  and  their  more  detailed  quantitative  study  has  necessitated 
the  developments  of  new  experimental  techniques  based  on  selective 
sampling  of  the  signal  during  bursts.  Although  different  criteria 
have  been  employed  by  different  investigators  as  to  what  should  signi¬ 
fy  the  occurance  of  an  "event",  all  are  in  substantial  agreement  as 
to  the  main  picture,  Willmarth  &  Lu  (1972)  discovered  that  60#  of 
the  Reynolds  stress  is  produced  when  the  velocity  in  the  sublayer  is 
less  than  the  mean.  Elackweid^r  &  Kaplan  (1972)  succeeded  in  deter¬ 
mining  the  average  instantaneous  velocity  profile  Just  preceding  burst 
for  different  locations  of  the  triggering  probe.  In  figure  9  is  shown 

the  instantaneous  profile  thus  obtained  before  burst  for  a  setting 

+  • 

of  the  triggering  probe  at  y  *15.  The  large  negative  velocity 
perturbation  up  to  20#  of  the  local  mean  and  the  strong  Inflexional 
character  is  apparent.  The  inflexional  character  of  the  velocity  profile 
preceding  burst  has  been  observed  in  all  investigations  of  the  phenom¬ 
enon  raid  suggests  a  strong  similarity  with  the  final  state  of  break¬ 
down  of  a  laminar  boundary  layer  to  turbulence  (Kiebanoff  et  al., 
1962),  Other  selective  sampling  data  have  been  presented  by  V/allace 
et  al.  (1972),  who,  extending  the  observations  by  Corino  and  Brodkey 
(1969),  found  that  the  ejection  and  sweep  stages  each  contributed 
about  seventy  per  cent  of  the  moan  stress  for  a  total  of  140  per  cent. 
The  excess  was  accounted  for  by  "interactions"  towards  and  away  from 
the  wall  each  producing  a  negative  turbulent  shear  stress  of  about 
twenty  per  cent  of  tho  mean.  In  between  ejections,  sweeps  and  inter¬ 
actions  the  flow  was  found  to  be  cf  essentially  laminar  character 
but  unsteady. 


Most  of  tho  experimental  data  on  the  turbulent  structure  of  a  nt”- 
Newtonian  fluid  reported  in  the  literature  so  far  have  pertained  to 
fluctuation  intensities.  From  these  it  is  immediately  apparent  that 
drag  reduction  is  not  accompanied  by  a  corresponding  decrease  in 
fluctuation  intensities.  In  fact,  the  measurements  reported  by  Virk 
et  al,  (1967)  showed  actually  an  increase  in  the  streamwise  fluctuation 
velocity,  Geyer  and  Metzner  (1969b)  on  the  ether  hand,  using  a  hydro¬ 
gen  bubble  technique,  found  a  decrease.  Development  of  las  ex'  anemo¬ 
meters  in  recent  years  has  alleviated  many  of  the  special  difficulties 
associated  with  measurements  in  polymer  solutions  employing  traditional 
techniques.  Measurements  by  Goldstein  et  al.  (1969)  gave  very  little 
difference  in  the  u-component  due  to  the  presence  of  on  additive. 

Rudd  0972)  found  a  noticeable  increase  in  the  streamwise  fluctuating 
component  near  the  wall  but  hardly  any  change  near  the  center  of  tho 
pipe.  (See  figures  10  and  11  reproduced  from  Rudd,  1972.)  Tho  trans¬ 
verse  (w-) component  w&3  decreased  near  the  wall- but  increased  above 
its  value  for  the  solvent  at  larger  distances. 

A  few  measurements  of  the  spatial  structure  of  ‘non-Newtonian  turbulences 
have  been  reported.  For  isotropic  grid  turbulence,  Fabula  (1966)  could 
not  detect  any  clear  effect  on  the  spectrum  from  polymer  addition, 
Friehe  {■:  Schwartz  (1970)  obtained  a  decrease  in  the  decay  rate  but 
no  ubstantial  change  in  the  spectrum.  Measurements  in  a  flat  plate 
boundary  layer  by  Kowalski  (1967)  showed  a  substantial  increase  in 
the  turbulence  microscale  near  tho  wall  with  as  much  as  by  a  factor 
of  two.  This  is  consistent  with  the  observations  from  mean  velocity 
measurements  that  tho  thickness  of  the  viscous  wall  layer  increases. 

On  basis  of  data  from  flow  in  rough  tubes  Virk  (1971a)  concludes  that 
the  maximum  increase  in  the  viscous  sublayer  due  to  polymer  additives 
is  about  2,5  in  qualitative  agreement  with  Kowalski's  result.  Results 


.27- 


from  correlation  measurements  have  only  recently  become  available. 
Fortune  &  Hanratty  (1972)  employing  an  electrochemical  technique 
determined  the  wall- stress  fluctuation  in  pipe  flow.  They  discovered 
a  strong  increase  in  the  transverse  (spanwise)  correlation  scale 
(see  figure  12)  as  well  as  a  much  slower  longitudinal  (streamwise) 
correlation  decay  (see  figure  13)  than  in  the  Newtonian  case.  A  re¬ 
duction  of  the  stress  fluctuation  amplitude  fox'  the  same  mean  wall 
stress  was  also  observed)  for  the  transverse  component  up  to  a  factor 
of  about  two,  but  less  for  the  axial  component.  Spectra  were  also 
determined  but  showed  no  dramatic  differences  from  their  Newtonian 
counterparts  when  normalized  with  wall  parameters.  However,  the  low- 
frequency  part  increased  somewhat  and  the  high  frequency  one  decreased 
with  increasing  drag  reduction. 

In  summary,  the  qualitative  picture  that  emerges  from  such  experimental 
results  is  that  polymer  addition  causes  an  increase  in  the  scale  of 
the  smallest  turbulents  eddies  near  the  wall  accompanied  by  a  reduction 
of  Reynolds  stress  and  an  associated  lowering  of  the  wall  friction. 
Turbulence  is  not  repressed  in  the  wall  layer  but  actually  enhanced. 

The  most  spectacular  change  in  the  turbulence  structure  is  in  the 
transverse  scale,  v/liich  can  Increase  with  up  to  a  factor  of  five 
(Fortuna  and  Hanratty,  1972),  The  reduction  of  streamwise  decay  of 
the  longitudinal  correlation  as  compared  to  the  solvent  points  to  on 
increase  in  the  life-time  for  the  big  eddies. 
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6,  Proposed  mechanisms 

A  number  of  suggestions,  all  of  a  qualitative  nature,  have  appeared 
in  the  literature  as  to  what  mechanism  might  be  the  primary  cause  of 
the  drag  reduction.  The  key  element  in  all  of  these  has  been  the 
thickening  of  the  Wall  layer,  first  observed  by  Elata  et  ei .  (1966). 
The  explanations  differ  widely,  however,  as  to  how  a  very  small  amount 
of  polymer  addition  could  interfere  with  the  turbulent  mixing  pro¬ 
cesses  so  as  to  produce  a  substantial  increase  in  the  scale  of  the 
smallest  stress-producing  eddies  near  the  wall.  There  is  even  differing 
opinions  on  such  a  fundamental  question  as  whether  an  explanation  on 
a  continuum  basis  is  sufficient.  Barenblatt  et  al.  (1965  and  1969) 
point  to  the  observed  appearance  of  large  agglomerations  of  molecules 
and  suggest  that  these  will  form  viscoelastic  particles  at  least  as 
large  as  the  Kolmogoroff  scale  of  turbulence  which  -  because  they  tend 
to  become  effectively  rigid  at  very  high  fluctuation  frequencies  -  will 
act  to  stabilize  the  small  viscous  eddies  in  the  wall  layer.  Their 
effect  on  the  inertial  large-scale  motion,  on  the  other  hand,  is  to  • 
lead  to  an  increased  intensity.  The  idea  that  entanglements  of  mole¬ 
cules  would  resist  the  motion  in  the  wall  laye*’  was  also  suggested  by 

Fabuia  ct  al.  (1966)  and  discussed  by  Lumley  (1967),  but  in  his  later 
review  article  (Lumley,  1969)  he  comes  cautiously  to  the. 

.conclusion  that  agglomerations  of  molecules  is  not  an  important  factor 

in  drag  reduction.  A  continuum  viewpoint  was  taken  by  Seyer  &  Met>;ner 

(1969a),  who  based  their  reasoning  on  the  model  of  large  stroamwlue 

eddies  for  the  wall  layer  proposed  by  Bakewell  &  Lumley  (1967)  and 

the  converted  Maxwell  constitutive  model  for  the  fluid.  According  to 

the  latter,  largo  values  of  normal  stresses  would  be  set  up  in  a  rapid 

elongaticnal  flow,  such  as  would  be  produced  when  the  two  Cake  well  - 

Lumley  counterrotating  vortices  force  the  fluid  up  from  the  wall 

between  them.  This  motion  would  that,  bo  strongly  resisted.  Pferninger 
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(1967),  following  Theodorscn's  (1955)  idea  that  the  turbulence  could 
be  modelled  as  a  distribution  of  horseshoe  vortices  formed  at  the 
wall  and  lifting  up  from  it  as  they  are  swept  do’/rstream,  proposed 
that  drag  reduction  was  associated  with  the  strong  resistance  to 
vortex  stretching  which  experiments  have  shown  to  be  exhibited  by 
some  (but  not  all)  drag  reducing  fluids  (see  also  Gadd,  1968).  The 
horseshoe  vortex  concept  lias  also  been  embraced  by  Black  (1969),  who 
however  thinks  that  stretching  of  the  more  or  less  longitudinal  vor¬ 
tices  in  the  horseshoe  "legs"  is  perhaps  not  as  important  as  the 
direct  stabilization  of  the  wall  layer.  Influence  on  the  energy  bud¬ 
get  of  the  turbulent  fluctuations  in  a  more  unspecified  manner  has 
also  been  invoked.  Astarita  (1965)  notes  tlvat  the  energy  dissipation 
ought  to  dccrcaso  in  a  viscoelastic  fluid  since  the  dynamically  pro¬ 
cess  will  become  more  elastic  an U  hence  more  reversible  at  high 
frequencies.  Thu3,  he  argue3,  the  drag  must  decrease  in  order  to 
preserve  the  energy  buler.ee.  Energy  arguments  were  also  used  by 
Gavins  (1964)  and  are  also  implicit  inVirk's  ( 19T1b>)  elastic  sub¬ 
layer  model.  Adsorption  of  a  resilient  layer  of  molecules  on  the  wall 

\ 

layer  has  also  been  proposed  as  contributing  to  the  thickened  wall 
layer  by  Bryson  et  al.  (1971),  but  no  explanation  has  been  suggested 
as  to  how  this  con  change  the  turbulent  processes. 


7»  Application  of  a  new  model 

To  make  any  quantitative  or  qualitative  assessments  of  such  proposed 
explanations  of  drag  reduction  is  presently  not  possible  because  of 
the  lack  of  an  adequate  theoretical  model  for  the  dynamics  of  shear 
flow  turbulence.  In  view  of  tho  extreme  complexity,  randomness  and 
nonlinearity  of  turbulence,  the  prospects  for  the  discovery  in  the 
near  future  of  a  reasonably  complete  such  model  arc  very  poor, 
indeed.  However,  in  view  of  the  very  drastic  effects  produced  by 
additives  one  could  reasonably  hope  tint  even  a  fairly  crude  theore¬ 
tical  model  -  if  retaining  the  essential  dynamics  of  the  turbulent 
processes  -  should  be  able  to  indicate  the  possibility  of  modifi- 
Cations  of  the  tui’buient  field  leading  to  drag  reduction,  at  least 
qualitatively.  We  shall  here  present  the  first  attempts  at  an  ex¬ 
tremely  simplified  such  model  which,  hopefully,  contain,  enough 
of  the  essential  dynamical  elements  of  the  turbulence  structure 
to  make  an  exploratory  study  of  the  effect  of  additives  possible. 


The  new  model  represents  an  outgrowth  from  ideas  in  classical  hydro- 
dynamic  stability  applied  to  a  fully  turbulent  flew,  but  it  also 
retains  Townsend's  (1956)  basic  concept  of  a  tv;o-laycr  structure. 

Tho  scale  of  wall  turbulence  being  roughly  proportional  to  distance 
from  the  wall,  one  would  expect  that  tho  significant  stress  and 
energy  producing  eddies  in  the  inner  "active"  layer  close  to  tho 
’/all  are  of  much  smaller  scale  than  those  in  the  "passive"  outer 
layer.  It  may  therefore  bo  meaningful  to  think  of  the  turbulence 
conceptually  as  consisting  of  a  coupled  motion  at  two  disparate 
scales:  a  primary  one  of  large  scale  (of,  say,  the  order  of  the 
boundary  l.r/er  thickness)  and  a  much  smaller  secondary  one  (of  the 


order  of  the  viscous  length  scale  v/u  ).  The  random  fluctuating 
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field  is  superimposed  on  a  mean  velocity  field,  which  to  a  good 
approximation  can  be  regarded  as  a  parallel  shear  flow.  Accordingly, 
if  the  fluctuation  velocities  are  small  compared  to  the  reference 
free  stream  speed  so  that  a  linearization  is  allowable,  a  simple 
model  equation  for  the  fluctuating  velocity  field  (when  Fourier 
transformed  in  x,  z  and  t)  would  be  of  the  Orr-Sommerfeld  type, 
and  a  decomposition  of  the  general  unsteady  motion  in  terms  of  the 
eigcnmodcs  of  the  Orr-Sommerfeld  problem  should  be  the  obvious 
way  to  proceed.  However,  nonlinearity  would  modify  the  picture  in 
three  important  ways. 

First,  the  primary  wave  field  (belonging  to  the  large-scale  motion) 
would  propagate  through  a  flow  filled  with  small-scale  eddies, 
hence  these  will  act  effectively  as  scattering  sources  for  the  large- 
scale  motion  producing  a  gross  effect  perhaps  much  like  that  of  an 
eddy  viscosity  (see  Reynolds  &  Hussain,  1y?2,  for  efforts  to  model 
this  effect).  Wc  shall  hero  ignore  this  effect  completely  finding 
some  Justification  in  doing  so  in  the  calculation  of  Landahl  (1967.) 
which  gave  propagation  and  decay  properties  of  the  large-scale 
pressure-producing  eddies  in  good  qualitative  agreement  with  ob¬ 
servations  without  including  eddy  viscosity. 


Secondly,  nonlinearity  is 
since  those  arc  damped,  ns 
obvious  excitation  sources 


required  to  excite  the  primary  waves 
shown  by  Landahl  (1967),  and  the  most, 
that  have  come  to  light  so  f pr  are 


the  turbulent  bursts.  As  the  recent  experiments  reviewed  above 
indicate  that  these  are  highly  localized  in  time  end  space,  the 
excitation  may  bo  modelled  ns  effectively  point  sources  on  the  scale 
of  the  primary  motion. 


Thirdly,  tho  small-scale  secondary  motion  sees  not  only  the  parallel 
mean  velocity  field  but  also  the  primary  fluctuating  field,  so  that 
the  secondary  waves  may  be  thought  of  as  riding  on  the  primary  waves 
in  very  much  the  same  way  as  capillary  water  w?  ves  on  a  large-scale 
gravity  wave.  On  the  scale  of  the  secondary  motion  the  mean  plus  the 
primary  field  will  be  a  slowly  varying  one  in  space  and  time  and  may 
therefore  locally  be  approximated  as  a  parallel  one. 


For  both  the  primary  and  the  secondary  motion  we  will  thus  have 
as  a  model  equation  one  of  the  Orr-Soinmerfeld  type  with  the  mean 
velocity  replaced  by  mean  plus  primary  in  the  latter.  The  secondary 
wave  could  thus  be  locally,  unstable,  for  example  because  of  in¬ 
flexion.  Now,  it  has  been  shown  recently  with  the  aid  of  kinematio 
wave  theory  (Landahl,  1972)  that  there  is  a  critical  condition  for 
a  locally  selfcxcited  secondary  wave  given  by 

cg  -  c0  (22) 

''here  c  is  the  group  velocity  of  the  secondary  wave  and  c  the 
«  0 

phase  velocity  of  the  primary  wave  on  which  it  is  riding,,  at  which 
the  secondary  wave  will  grow  to  very  large  amplitudes,  Tho  reason 
for  this  is  that  tho  group  of  secondary  waves  at  this  condition 
experiences  space-time  focusing  (i.e,,  neighboring  groups  of  waves 
catch  up  with  each  others  ana  hcnco  "compress"),  and  that  the  wave 
group  becomes  trapped  in  a  region  of  the  primary  wave  (generally 
near  its  crest),  where  it  can  continue  to  amplify  (see  figure  It) . 

It  is  therefore  hypothesized  that  (22)  should  identify  tho  condition 
at  which  the  breakdown  of  the  primary  waves  into  strong  secondary 
small-scale  notion  would  take  place.  Application  to  the  data  from 
the  transition  experiments  by  Ii.lcbr.noff  et  al.  ( I962)  demonstrates 
that  the  condition  is  very  closely  satisfied  at  the  appearance  of 


the  fiiot  "spikes"  in  the  final 


r.-tfgc  of  transition  of  a  laminar 
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boundary  layer,  and  that  it  gives  good  agreement  with  the  measured 
frequency  of  secondary  motion.  In  view  of  the  observed  similarities 
between  breakdown  in  the  final  stage  of  transition  and  turbulent 
bursts  (as  pointed  out,  for  example,  by  l/illmarth  and  Lu,  1972)  we 
shall  hypothesize  that  the  bursts  are  the  results  of  tho  attaiatment 
of  local  criticality  in  tho  primary  field,  and  that  the  subsequent 
violent  motion  will  cause  excitation  of  new  primary  waves,  (The 
model  is  illustrated  conceptually  in  figure  15.)  Admittedly,  the 
model  is  highly  speculative  at  this  stage  and  contains  many  diffi¬ 
cult  details  which  are  not  yet  worked  out  and  which  will  require 
a  great  deal  more  study,  such  as  how  randomness  could  be  incorpo¬ 
rated  into  The  model,  and  how  primary  waves,  which  probably  are 
highly  oblique,  could  interact  so  as  to  produce  a  critical  condition. 
Nevertheless,  fer  the  purpose  of  studying  the  effects  of  additives, 
a  great  deal  could  perhaps  be  learned  from  the  propagation  charac¬ 
teristics,  in  particular  the  stability  characteristics,  primary 
and  secondary  motion,  Accordingly,  some  numerical  stability  calcu¬ 
lations  have  been  carried  out  for  velocity  profiles  that  may  bo 
typical  of  tho  mean  and  instantaneous  primary  motion  for  both  a 
Newtonian  situation  and  one  with  drag  reduction.  These  calculations 
were  performed  by  P.  Dark  of  the  Royal  Institute  of  Technology, 
Stockholm,  and  a  detailed  account  of  them  will  be  given  later. 

The  three  rheological  models  that  were  used  were  the  ones  described 
above;  tho  ordinary  Newtonian  one,  the  convected  Maxwell  one  and 
the  Batchelor  -  Hinoh  model  for  a  suspension  of  microscopical  rigid 
rods.  Three  different  velocity  profiles  were  investigated,  two 
qualitatively  representative  of  tho  mean  motion  for  a  Newtonian 
and  a  drag- reducing  fluid  in  a  boundary  layer,  aid  one  broadly 
similar  to  the  inf  J.  ox  tonal  profile  appearing  in  the  sublayer  .Jur-t 
before  bre'\kdo..n.  All.  throe  arc  illustrated  in  figure  16,  'J’hc  mean 
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profllo  for  tho  Newtonian  caoo  wa3  constructed  from  Reiohhardt' s 
representation  for  the  velocity  distribution  in  the  constant- stress 
layer  coupled  with  Coles'  wake  law  profile  for  the  outer  portion# 

The  drag- reduction  velocity  profile  was  obtained  in  essentially 
the  same  'way,  except  that  tho  viscous  wall  layer  thickness  para¬ 
meter  in  Rcichhardt's  expression  wa3  increased  to  correspond  to  the 
maximum  wall  layer  thickness  possible  according  to  the  Seyer  - 
Metzner  (1969b)  correlation.  The  same  plate-length  Reynolds  number 
was  employed  as  for  the  Newtonian  profile.  By  comparisons  of  re¬ 
sults  for  the  two  profiles  with  tho  same  constitutive  model  one  will 
thereby  be  able  to  isolate  the  effect  on  the  stability  characteris-  * 
tics  of  the  change  of  profile  from  an  ordinary  turbulent  one  to 
a  more  laminar-looking  one. 

The  velocity  profile  chosen  to  modol  the  instantaneous  wall  layer 
breakdown  profile  was  of  the  simple  tanhyp-variety  with  tho  distance 
of  tho  inflexion  point  to  the  wall  adjusted  30  as  to  be  roughly  in 
the  same  ratio  to  the  internal  shear  layer  as  in  the  experimental  . 
profile  determined  by  Blackwclder  and  Kaplan  (1?Y2).  Tills  vias  selec¬ 
ted  in  preference  to  a  direct  numerical  fit  to  the  measured  profile 
In  order  to  avoid  undue  numerical  difficulties  that  would  arise 
because  of  the  fourth  derivative  of  U  appearing  in  the  equation 
for  tho  connected  Maxwell  model.  Such  simplifications  should  not 
be  serious  at  the  present  stage,  however,  where  only  qualitative 
indications  of  a  strong  effect  arc  being  sought. 

Results  from  tho  primary  wave  calculations  are  shown  in  figure  17. 
Following  Landahl  (I967)  these  were  obtained  in  spatial  form  with  tho 
complex  wave  number  a  =  <yu  +  io'  determined  as  function  of  the  real 

frequency  w  (made  dimension]. ess  with  the  "outer"  reference  quanti- 
i. inf;  a  nml  11. 1.  'tho  imcprinarv  t>art.  fives  a  measure  of  tho 
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rate  of  decay  with  streamwise  distance  for  a  wave  of  given  frequency. 
The  results  obtained  with  the  Maxwell  fluid  model  show  that  for  the 
lower  frequency  range,  ui  <  1.5,  a  moderate  fluid  elasticity  (non- 
dimensional  relaxation  time  0  less  than  unity)  causes  very  little 
change  in  relative  to  the  Newtonian  model  for  the  same  velocity 
profile.  There  is  an  increase  of  orr  relative  to  the  value  for  the 
Newtonian  profile  which  is  entirely  a  profile  effect.  At  higher 
frequencies  a  destabilizing  effect  is  seen  which  at  w  »  3  is 
strong  enough  to  lead  to  a  complete  loss  of  stability.  Calculations 
fur  larger  values  of  0  (not  shown  here)  give  frequency  regions 
in  which  qtj.  becomes  negative,  i.e,,  the  wave  is  unstable.  It  is 
difficult  to  Judge  the  significance  of  such  results  as  the  model 
employed  requires  that  all  primary  waves  are  stable,  otherwise  the 
assumption  of  statistical  stationarity  has  to  be  given  up  (Lindahl, 
1967). 

With  the  Batchelor  -  Hinch  constitutive  model,  on  the  other  hand, 
the  decay  rate  o^  is  increased  over  the  value  for  the  Newtonian 
case  by  about  a  factor  of  2, 5.  At  a  given  non-dimensional  frequency 
based  on  outer  variables  one  would  thus  expect  to  see  a  correspon¬ 
dingly  more  rapid  decay  of  the  cross-power  spectrum  of  the  turbulent 
fluctuations  with  streamwiso  distance.  Such  a  result  might  seem  to 
bo  at  variance  with  the  experiments  of  Fortuna  4  Hanratty  (1972), 
which  showed  a  much  reduced  decay  of  the  streamwise  correlation 
of  the  wall  stress  fluctuations  with  stream',.'! se  distance-  as  compared 
to  the  Newtonian  case  for  tho  came  mean  wall  stress.  However,  sinco 
the  correlation  is  the  inverse  transform  of  the  cross-power  spectrum, 
the  reduction  in  the  correlation  decay  rate  may  simply  be  due  to  a 
shifting  of  the  whole  spectrum  towards  lower  frequencies,  as  the 
experiments  Indeed  indicate  will  take  place  during  drag  reduction. 
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Results  of  the  stability  calculations  for  the  secondary  disturbance 
with  the  model  inflexional  profile  selected  are  presented  in  figure 
18.  In  these  calculations  the  fairly  high  Reynolds  number  of  Re  ■  ^600 
was  chosen  in  order  to  make  the  Newtonian  stability  characteristics 
become  essentially  the  inviscid  ones.  For  a  study  of  the  shortest 
scale  in  a  Newtonian  wall  layer,  perhaps  a  more  realistic  value  of 
the  Reynolds  number  would  have  been  less  than  1000;  however,  in 
order  to  bring  out  more  clearly  any  stabilization  potentialities 
of  the  other  rheological  models,  a  high  Reynolds  number  was  chosen. 
From  the  results  presented  in  figure  17  it  becomes  clear  that  a 
moderately  iiigh  elasticity  of  the  fluid  (0  »  0(1))  has  only  a  small 
destabilizing  influence.  However,  with  the  Batchelor  -  Hinch  model 
a  strong  stabilizing  effect  is  demonstrated  including  the  possi¬ 
bility  of  a  complete  stabilization  of  the  inflexional  profile  for 
sufficiently  high  values  of  the  parameter  B,  In  particular,  it  is 
observed  that  the  strongest  stabilizing  effect  is  found  for  the 
largest  wave  numbers,  so  that  the  smallest  scale  of  unstable  motion 
is  increased.  Thi3  would  he  expected  to  lead  to  an  increase  in  the 
size  of  the  smallest  stress  producing  eddies  and  a  corresponding 
increase  in  the  wall  layer  thickness.  The  values  of  5  used  in  the  ' 
calculations  may  perhaps  seem  Iiigh,  but  they  are  not  unrealistic, 
as  they  give  a  measure  of  the  ratio  of  clongational  to  shear  visco¬ 
sity  of  the  same  order  as  found  in  the  experiments  by  Metznor  & 

Metzncr  (1970)  •  Furthermore,  it  is  seen  that  tho  parameter  enters 
into  the  basic  stability  equation  (22)  in  the  ratio  B/Re,  so  that 
for  a  lower  Reynolds  number  the  value  of  B  required  for  stabilization 
becomes  proportionally  lower. 
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The  simplified  turbulence  model  presented  thus  shows  that  fluid 
elasticity  as  modelled  by  the  Maxwell  constitutive  equation  causes 
destabilization  of  both  primary  and  secondary  motion.  The  Batchelor 
Hinch  constitutive  model,  on  the  other  hand,  is  found  to  produce 
a  strong  stabilization  of  both  primary  and  secondary  motion.  The 
stabilizing  effect  on  the  secondary  motion  is  probably  of  most 
significance  for  the  drag  reduction  phenomenon  since  the  associ¬ 
ated  increase  in  scale  of  the  unstable  secondary  disturbances  could 
be  expected  to  lead  to  a  thickening  of  the  wall  layer. 


8,  Concluding  remarks 


The  first  basic  and  as  yet  unresolved  question  in  the  puzzling  phenom¬ 
enon  of  drag  reduction  due  to  polymer  addition  is  whether  a  satis¬ 
factory  explanation  based  on  continuum  mechanics  is  possible  or 
whether  one  instead  would  have  to  invoke  processes  on  the  non-continuum 
macromol ccuiar  level  such  as  proposed  by  Ba^enblatt  et  al.  (1965  &  1969). 
From  a  heuristic  point  of  view,  a  theory  based  on  non- continuum  con¬ 
cepts  will  be  quite  difficult  to  quantify  and  the  theoretician  would 
therefore  naturally  be  quite  reluctant  to  part  with  the  continuum 
approach.  Although  the  Initial  efforts  towards  a  simplified  mechanistic 
theory  of  wall  turbulence  presented  here  have  not  yet  led  to  definite 
conclusions  in  this  regard,  they  have  at  least  pointed  to  a  possible 
continuum  mechanical  effect  due  to  the  polymer  additive  which  could 
be  highly  effective  even  when  the  molecular  size  is  much  smaller  than 
the  scale  of  fluctuating  motion.  The  large  frictional  resistance  to 
streamwise  straining  perturbations  caused  by  elongated  particles 
aligned  in  the  mean  flow  direction  was  shown  to  have  a  strong  stabili¬ 
zing  effect  on  inflexional  Instability,  Since  there  is  both  experi¬ 
mental  and  theoretical  evidence  that  such  instability  is  involved  in 
the  creation  of  the  stress  producing  turbulent  bursts  observed  to 
occur  intermittently  in  the  wall  layer,  such  stabilization  could  be 
expected  to  lead  to  a  reduction  of  the  turbulent  stresses  near  the 
wall  and  an  associated  thickening  of  the  wall  layer.  Extended  polymer 
macromol ocules  mostly  aligned  in  the  mean  flow  direction  of  the  tur¬ 
bulent  velocity  field  could  bo  expected  to  act  qualitatively  like  the 
rigid  rods  employed  in  the  simplified  fluid  model,  at  least  during 
the  extension'll  phase  of  the  secondary  fluctuation  cycle.  Tint  the 
polymer  molecules  could,  have  a  rhcolc-gtcal  effect  of  tills  kind  is 
indicted  inc.t--v;i-jy  by  the  very  hi'h  •/clues  of  the  ratio  of  donga- 
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tional  to  shear  viscosity,  of  the  order  1o\  measured  by  Metzner  ft 
Metzner  (1970).  Such  high  values  are  consistent  with  what  Batchelor's 
^971 )  simplified  theory  would  give  for  rod.'-:  having  aspect  ratios  like 
those  of  the  extended  polymer  molecules  of  typical  drag  reduction 
additives  in  dilute  solutions. 


As  yet  no  direct  experimental  demonstration  of  the  presence  of  fully 
extended  macromolocu]  ar  threads  during  drag  reduction  has  been  given. 
Complete  extension  would  not  bo  attainable  in  a  laminar  shear  .  .  '  .  • 
fjow,  but  the  fluctuating  irrotnlional  strains  present  in  the  turbulent 
field  could  possibly  cause  sc..ie  ''.oils  to  stay  extended  a  large  fraction 
of  time.  Intuitively,  or.o  would  feel  that  the  onset  phenomenon  is  a 
manifestation  of  the  extension  of  the  macrorrol ocular  coils,  although 
the  observation  (Virk  et  al.  ,  1967)  that  onset  data  best  correlate 
with  a  length  hypothesis  based  on  the  ratio  between  coil  radius  and 
the  smallest  scale  of  turbulence  has  not  yet  found  a  satisfactory 
theoretical  explanation. 


Reproduced  from 
best  available  copy. 


The  calculations  based  on  the  Maxwell  fluid  model  give  quite  clear 
indications  of  the  effects  of  fluid  elasticity  on  fluctuation  field 
stability,,  J5oth  the  large-scale  primary  motion  and  the  small-scale 
secondary  one  were  found  to  bn  destabilized  by  moderate  amounts  of 
elasticity.  V.'Ith  the  model  cr .ployed  there  to  thus  no  indication  that 
viscoelasticity  as  much  could  ho  the  caus*.  of  drag  reduction. 


The  explanation  that  drag  reduction  is  associated  wi  th  the  high  clonga- 
tional  viscosity  produced  by  aliened  uncoiled  macromolecular  threads 
is  broadly  in  accord  with  Virk' s  0971b)  findings  that  the  drag  curve 
slope*  increment  6  is  approximately  proportional  to  the  square  root 
of  concentration  and  the  three  half  power  of  the  number  of  backbono 


parameter  D  is  proportional  to  nl3;  (ignoring  the  logarithmic  term) 
where  n  is  the  number  density  and  1  the  length  of  the  aligned  particles. 
The  experimentally  determined  variation  with  concentration,  i.e.,  of  n, 
to  the  one  half  power  would  thus  indicate  the  effect  to  be  proportional 
to  (nl3)  /  ,  which  thus  yields  a  variation  with  the  bond  length  b 
and  number  of  backbone  links  ::  as  b3/2  N3/2.  Thus,  proportionality 
to  c1/2  would  bo  consistent  with  the  H3/2  variation*  The  two  models 
differ  as  to  the  variation  with  molecular  chain  bond  length,  however. 

On  the  basic  of  an  estimate  of  the  turbulent  strain  energy  stored  in 
the  clastic  molecules,  Virk  (1971b)  predicted  a  proportionality  to  b3, 
in  contrast  to  the  b^/2  dependence  inherent  in  the  (nl3) 1/2  variation. 
The  nigh  drag  reduction  values  exhibited  by  guar  gum  are  more  in  line 
with  the  b3  variation  propose,!  by  Virk  than  with  the  proportionality 
to  b  tending  to  givo  .  uroort  for  Virk's  hypothesis  that  elasticity 
is  an.  important  molecular  property  in  drag  reduction.  In  the  Batchelor  - 
ILnCii  model  only  he  t.-bal  length  of  the  particle  plays  a  role;  its 
hydrodynamic  effect  c  '  ndo  over  a  sphere  of  diameter  1  and  the  quan¬ 
tity  nl3  thus  measures  the  number  of  such  influence  spheres  per  unit 
volume.  The  high  values  of  B  that  would  be  typical  for  dilute  drag 
reducing  fluids  -  of  the  order  of  103  if  the  macromolccules  were  fully 
extended  -  indicate  that  tho  hydrodynamic  interaction  between  the  ex¬ 
tended  macromolecular  threads  should  be  very  strong  (i.e.,  their 
spheres  of  hydrodynamic  influence  would  overlap  strongly).  In  order 
for  tho  solution  to  bo  dilute  in  the  rheological  sense  s0  that  the 
molecules  take  up  only  a  small  fraction  of  the  available  volume,  they 

have  to  have  extremely  high  length  to  diameter  ratios  of  tho  order 
103  or  moro. 


Tills  was  pointed  out  to  me  by  E.J,  Hinch. 
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However,  there  is  recent  experimental  evidence  that  it  may  not  only  be 
the  length  of  the  extended  molecule  that  ia  of  importance  ' 
for  drag  reduction.  Hand  &  Williams  (1971)  measured  the  drag  reduction 
effectiveness  of  polyacrylic  acid  and  of  calf  thymus  DNA  at  various 
pH  levels  of  the  solution.  They  found  that  the  effectiveness  was  con¬ 
siderably  larger  for  pH  ranges  in  which  the  molecules  are  believed 
to  have  a  fatrxy  compact  helical  structure.  For  pH  ranges  where  the 
molecules  are  supposed  to  be  extended  and  behave  essentially  as  rigid 
rods,  the  drag  reduction  was  found  to  be  considerably  lower.  Results 
such  as  these  .re  very  difficult  to  explain  on  basis  of  any  continuum  * 

modelj  but  possibly  entanglements  and  agglomerations  could  be  an 
important  factor. 

It  should  be  pointed  out  that  stabilization  of  the  instantaneous 
inflexional  velocity  profile  before  burst  -  as  indicated  by  the  cal¬ 
culations  for  a  suopension  of  rod-like  particles  -  is  only  one  of 
several  conceivable  mechanisms  that  could  produce  a  decrease  in  the 
rate  or  intensity  of  tne  bursts  and  thereby  cause  a  lowering  of  tur¬ 
bulent  stress.  An  alteration  of  propagation  characteristics  of  the 
secondary  wave  group  in  relation  to  the  primary  disturbance  field 
could  presumably  have  ouch  an  effect  according  to  the  recently  pro¬ 
posed  model  for  breakdown  (Londahl,  1972),  but  the  model  has  not  yet 
been  sufficiently  developed  to  allow  any  conclusions  to  be  drawn  re- 
gardiiig  the  possibility  of  an  effect  of  this  kind.  Conceivably,  the 
enhunccd  clongational  viscosity  could  also  have  a  strong  influence 
on  the  structure  of  the  longitudinal  vortices  appearing  in  the  wall 
layer,  in  particular,  a  strong  resistance  to  vortex  stretching  would 
bo  expected.  However,  a  quantitative  ^del  to  ’  investigate  this  effect 
is  likely  to  be  quite  difficult  to  set  up. 
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As  pointed  out  in  the  Introduction,  perhaps  the  most  interesting 
aspect  of  the  research  on  the  drag  reduction  effect  of  polymer  addi¬ 
tives  is  vrtiat  new  knowledge  can  be  gained  from  it  regarding  the  nature 
of  turbulent;.  That  one  now  has  a  very  effective  tool  with  which  the 
turbulent  processes  can  be  manipulated  should  open  up  new  possibilities 
for  further  experimental  and  theoretical  research  on  the  mechanics  of 
turbulence.  For  example,  such  experiments  as  those  by  Wells  &  Spangler 
(1967)  on  the  effectiveness  of  injection  of  polymer  at  the  wall  demon¬ 
strate  quite  unequivocally  that  wall  shear  flow  turbulence  is  control¬ 
led  by  processes  in  a  layer  very  close  to  the  wall. 


An  important  engineering  motivation  for  turbulence  research  is  to 
find  new  and  effective  methods  for  control  of  turbulence.  From  the 
research  on  polymer  additives  reviewed  here  it  becomes  apparent  that 
to  Influence  skin  friction  it  is  the  very  small  -scale  stress  producing 
eddies  near  the  wall  which  may  need  to  be  interfered  with  directly. 
This  may  put  very  difficult  constraints  on  the  dynamical  response 
characteristics  of  any  mechanical  system  employed  for  drag  modifi¬ 
cation,  but  the  search  for  efficient  such  systems  should  bo  an  ex¬ 
tremely  interesting  challenge  from  both  a  scientific  and  engineering 
point  of  view. 


<xi  mu..,  jij.„:>i.od  to  j'-.J.  Iftncn  of  Cambridge  University  fer  many 
interesting  and  useful,  discussions  on  the  subject  of  drag  reduction 
which,  as  is  evident  from  the  paper,  spawned  many  of  the  new  ideas 
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of  .'.usjK'nsin  i  .. 
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ADDITIVE -FLUID  COMBINATION'S  SHOWING 


DRAG  REDUCTION 

ADDITIVE  FLUID 

Polyethyleneoxide  (PEO,  Polyox  )  water 
Polyacrylamide  (PAM,  Separcn) 

Guar  gum  (GGM) 


Polymethylmethacrylate  ( PMM )  monochlorobenzene 
. .  (Toms  ,19^8) 

Vaf  ious  soaps  water,  toluene 

Nylon  fibres  water 

( Bobkowicz  $  Gouvin,  1365  ) 


licure  1,  uouic  arloi  Live- fluid  .syutoiriO  exhibiting  drag  induct'*  or. 


i 


Ficure  2.  Pipe  flew  drag  data  for  solutions  of  oolyethylencojd.de 
(from  Paterson  &  Abernathy,  1970) 


Fitiuro  3.  Vhe  i;i?.yitnur.i  drag  reduction  nc.yr;>tote 
(from  V.lr’c  at  al.,  1970) 


•Figure  4.  Expc.rimi.ntal  mean  velocity  profiles  ourinc  dm?  reduction, 
(from  Virk,  I'iyibl .  A  A  t  Plata  cc  al.  (1966); 
ti  in  Goron  &  Horbury  0967);  V  Patterson  fc  Florez  (1969; 

S>  O  Ecyar  &  Mctzr.er  (1969b)  j  o  O  VI  rk  ct  al,  (I96Y) 


J'Mf.ni'e  9.  f/>r;.avl  l.hmlc  .Intercept  function  for  dm;;  reduction 
velocity  profile;;  bey  o  r  <>;  Me  tznov,  1969b) 

Data  obtained  with  nclvaervlrrildo  in  water. 
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FiGurc  6.  Intrinsic  slope  Increment  factor  jr  rs  function  of  number 
of  backbone  links,  H  (from  Vivk,  1971b)  O  guar  gum,  water; 
0* polyacrylamide,  water; polycieisoprone,  toluene; 

Pa  polye’-ime thylslioxan e,  toluene;  OOpplyethyieneoxldc, 
water,  benzene  0.6  m  A  polylsobutylene,  eyelehoxsne; 

v  po lym etliy  1  tn o thac l y  1  at o ,  chlorobenzene. 


Figure  7,  JVvtohclor's  (1Q71 )  model  cf  a  suspension  of  aligned 
rigid  rods. 


Figure  9.  Instantaneous  velocity  distribution  Just  preceding 
burst  Measured  by  Bln.ckwelder  &  Kaplan  (1972) 
Trigger  sot  at  y+  =  15t 
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F.lfsurc  12.  Spatial  correlation  coefficient  of  atreamwise  fluctuating 
v/all  stresses  (from  For  tuna  &  Hanratty,  1977) 

water;  -m~o-  solution  of  polyacrylarnld  giving  65# 
drag  reduction. 
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W'.-m'c:  O.  Spatial  correlation  coefficient  of  transverse  fluctuating 
wall,  stresses  (fnrn  For  Luna  .1 11, -n  rally,  1972).- —  water; 
0  solution  of  polyacryltaJ.do  giving  to. dra «  ur.lmvM™ 


SECONDARY 
(Small  scale  ) 
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